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Abstract 

Some known constraints on Renormalization Group flow take the form of inequalities: in even 
dimensions they refer to the coefficient a of the Weyl anomaly, while in odd dimensions to the 
sphere free energy F. In recent work |lj it was suggested that the a- and F-theorems may be 
viewed as special cases of a Generalized F-Theorem valid in continuous dimension. This conjecture 
states that, for any RG flow from one conformal fixed point to another, Fuv > Fir, where F = 
sin(7rd/2) log Z S d. Here we provide additional evidence in favor of the Generalized F-Theorem. We 
show that it holds in conformal perturbation theory, i.e. for RG flows produced by weakly relevant 
operators. We also study a specific example of the Wilson-Fisher 0(N ) model and define this CFT 
on the sphere paying careful attention to the beta functions for the coefficients of curvature 

terms. This allows us to develop the e expansion of F up to order e' 5 . Pade extrapolation of this 
series to d = 3 gives results that are around 2 — 3% below the free field values for small N. We also 
study RG flows which include an anisotropic perturbation breaking the 0(N ) symmetry; we again 
find that the results are consistent with Fuv > Fir. 
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1 Introduction and Summary 

A well-known set of constraints on Renormalization Group in d-dimensional relativistic Quantum 
Field Theory (QFT) takes the form of inequalities: if an RG trajectory leads from a short-distance 
unitary Conformal Field Theory (CFT) to a long-distance one, then a certain positive quantity 
defined on the space of CFTs decreases. In even dimensions d, using the modern terminology such 
theorems may be called the a-theorems. The Weyl anomaly equation takes the form 


{Tp^{-l) d / 2 a E d + Y J C i Ii, 

i 


(l.i) 


where Ej is the Euler density term, which is present in all even d > 2, and Cj are the coefficients of 
other Weyl invariant curvature terms. The universal Weyl anomaly coefficient a may be extracted 
from the free energy on a sphere of radius R: F = — log Z S d = (—l) rf//2 a log R. In even d, the RG 
inequalities take the form 12—[7] 

auv > am • (1.2) 


In d = 2 this is well-known as the c-theorem |2 , because there is only one Weyl anomaly coefficient, 
and it is standard to define c = 3a. In d = 4 a non-perturbative proof of (1.2) was found in 
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long after the early work of |3|8]. Very recently [6||7], there was major progress towards establishing 
(1.2) for supersymmetric flows in d = 6, building on the earlier work including [9,10 . 


In odd dimensions the situation is quite different because there are no Weyl anomalies. This 
actually has some advantages: for a CFT, the sphere free energy F = — log Z S d is a finite, radius 
independent quantity, which has no ambiguities because there are no Weyl invariant terms con¬ 
structed purely out of the curvature tensor. In d = 3 it was conjectured 111 that the RG inequality 
takes the form Fry > Fir. This F-theorem can be equivalently formulated in terms of the entan¬ 


glement entropy across a circle 12,13 . A proof of the three-dimensional F-theorem has been found 


using properties of the entanglement entropy in relativistic theories [14] (see also [15]). For other 
odd dimensions, it was conjectured 16 that the RG inequality takes the form Fry > Fir, where 


F = (—l)( d 1 ^/ 2 log Z S d J^] Supporting evidence for this conjecture included a calculation of the 
change in F for the RG flow produced by a weakly relevant operator of dimension d — e 


16 


Such 


calculations, as well as their analogues in even dimensions which apply to the change in a 
can be carried out perturbatively in e. 

The similarity between the RG inequalities in even and odd dimensions (both of them can be 
phrased in terms of the free energy on S d or, equivalently, in terms of the entanglement entropy 
S d ~ 2 12 ) has led to the idea that they are special cases of the RG inequality valid in 


across 


continuous dimension 111] 2 


Fry > Fir , F = sin(7rd/2) log Z S d = — sin(7rd/2)F . 


(1.3) 


In even dimensions, the factor sin(7rd/2) cancels the pole present in F, and we have F = 7ra/2; 


therefore (1.3) reduces to the a-theorem. In odd dimensions, this definition reduces to that proposed 
in 


16 . Thus, the “Generalized F-Theorem” (1.3) smoothly interpolates between the a-theorems 


in even d and the F-theorems in odd d. 

In |1| several pieces of evidence were provided in favor of the Generalized F-Theorem J^] In 
free conformal scalar and fermion theories, F is positive for all d. For example, for a conformally 
coupled scalar it is 


Fs(d) = 


1 


r (i + d) Jo 


1 (d \ (d 

du u sin mil — + u I [ ~ u 


(1.4) 


which is a smooth, monotonically decreasing function of d. Therefore, (1.3) is valid for the RG 


flows produced by the scalar or fermion mass terms - such flows lead to trivial theories where 
F = 0. Other tractable RG flows include those produced by double-trace operators in large N 


17 


Tn the special case d = 1 this conicides with the g-theorem for boundary conformal field theory 
“Introduction of the factor sin(7rd/2) is just a convenient choice; a further multiplication of F by a positive function 
of d does not change the conjectured inequality. 


3 Some holographic evidence for the Generalized F-Theorem was also provided in 18 . 
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theories. As is well-known 19,20 , turning on a relevant perturbation O?, which is the square of of 


a primary scalar operator of dimension A < d/2, makes a large N CFT flow to another CFT where 
the corresponding operator has dimension d — A + 0(1/N). This produces the following change in 


F0 


21 


Fir ~ F\jy = 


1 


r(i + d) Jo 


A 2 / d 

du u sin itu 1 1 - + u 


--u)+0{l/N). 


(1.5) 


This is negative in the entire range (d — 2)/2 < A < d/2 where the UV and IR CFTs are unitary, 
establishing consistency with the Generalized E-Theorem. However, when A is sufficiently far below 


the unitarity bound (d — 2)/2, then the inequality (1.3) is violated (for d = 3 this was discussed 


in [22]). This provides a simple explicit example of how (1.3) may be violated in non-unitary 
theories]!] 

Explicit applications of the a-theorem in d > 2 have been mostly to supersymmetric RG flows 
(see, for instance, §0 24/28]), because few interacting non-supersymmetric CFTs are known in 
d = 4 and none in d = 6. On the other hand, in 2 n — e dimensions there is a multitude of non- 


supersymmetric RG flows connecting perturbative fixed points [29-35 . This opens the possibility of 


many tests of the Generalized E-Theorem in 2 n — e dimensions, and some of them were carried out 
in |T]. In this paper we extend these calculations, focusing mostly on the Wilson-Fisher fixed point 
in 4 —e dimensions 29 , which is the 0(N) symmetric theory of N real scalar fields </>*, i = 1,..., N, 


with interaction ^(0*0*) 2 . For large N , this is a double-trace operator with A = d— 2. Furthermore, 
the fact that the coupling constant at the IR fixed point is of order e allows one to develop the e 
expansions for the critical exponents |29|; this works well for all values of N including N = 1, i.e. 
for the Ising model 36 . Fortunately, E is also amenable to e expansion using renormalization of the 
0(N ) model on the sphere S 4_e . In [l the contributions of interactions to E were determined up to 
0(e 4 ). The leading interaction correction, which is of order e 3 , did not require renormalization and 
was quite straightforward. However, at the next order one needs to renormalize the theory including 
the effects of the terms quadratic in the curvature. In section [2] we elucidate the definition of the 
Wilson-Fisher fixed point on S 4 ~ e following [8 37-40]. We find that it is important to define the 
IR theory by setting all the beta functions to zero, including the beta functions for the coefficients 


In non-integer dimensions all theories appear to exhibit a subtle form of non-unitarity 


23 


but hopefully this 


does not invalidate the Generalized T-Theorem for theories that are dimensional continuations of unitary theories in 
integer d. 
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of the curvature terms. Applying this procedure up to order e 5 we find 


P _ A 7 tN(N + 2) 3 nN(N + 2)(131V 2 + 3701V + 1588) 4 

° (N) ^ s ^ ' 576(1V + 8) 2 6912(1V + 8) 4 

+ 4 J^y + 2 ) (l0368(iV + 8)(51V + 22)C(3) - 647A 4 - 321521V 3 - 6065761V 2 - 39395201V 

+30tt 2 (A + 8) 4 - 8451008) e 5 + C>(e 6 ). 

As a byproduct of our analysis of the 0(N ) model on S 4 ~ e , we extend the previous results [8 
and determine the beta function for the Euler density term to order A' 5 : 


( 1 . 6 ) 


38 


40 


Pb = eb- 


N 


N(N + 2)(3N + 14) 4 N(N + 2)(3(24 + 7N) + 4(51V + 22)C(3)) 5 


360(4vr) 


48(4tt) 


10 


-a 4 + ~ v ~ 1 ~~ /,v-// a 5 +0(a 6 ). 

(1.7) 

As far as we know, the 0( A 5 ) term has not appeared in the previous literature. We also present 
a calculation of the sphere free energy and curvature beta functions in the most general quartic 
scalar field theory in Appendix [Aj 


A constrained Pade approximation of the series (1.6) for N = 1, using the boundary condition 


that the 2-d Ising model has F = vr/12, corresponding to c = 1/2, gives that in 3-dimensions 
Ei S ing/Es ~ 0.976. This is consistent with the E-theorem, but the surprise is how close Ei s ; ng is 
to the free field valuej^] For comparison, we note that the scaling dimension of q ~ 0.5182 is 
around 3.6% above the free held value. The coefficient of the stress tensor 2-point function, ct, is 
also known to be close to the free held value for the 3-d Ising model: c^ Ismg /cy ~ 0.9466 


41 


42 


In section [4] we also provide a check of the Generalized E-Theorem for weakly relevant pertur¬ 
bations. When an RG how in d dimensions is sourced by an operator O of dimension d — e then 
there is a nearby IR hxed point, and the change in E is 


Fir - Euv = - 


C 3 


e 3 + 0 (e 4 ) 


( 1 . 8 ) 


r (1 + d) 3Cf 

where C 2 and C 3 are the coefficients of the two- and three-point functions of O in the UV CFT, 


(4.2). In odd d this result agrees with 16 , and in even d with the change in a-anomaly computed 


in j5j. For a unitary CFT, C 2 > 0 and C 3 is real. So we hnd that the Generalized E-Theorem holds 
to leading order in conformal perturbation theory for all d. We have also extended the conformal 
perturbation theory approach to determine the term of order e 4 in specihc models. For example, 


for the 0 (N ) model it reproduces this term in ( 1 . 6 ). 


In section |4~2] we generalize the conformal perturbation approach to the case of several weakly 


relevant operators. We obtain a concise formula (4.39) for the leading order change in E. As a 


s The numbers we get from the Pade approximants are close to the the estimates in 1 which were obtained without 
the use of resummation. 
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specific application, we study the O(N) model with an extra “cubic anisotropy” operator Yi^t 
which breaks the 0(N ) symmetry. The conformal perturbation approach allows us to calculate the 
F to order e 3 at the different fixed points of this theory with two coupling constants. The results 
are found to be in agreement with the direct approach of section [2] which uses renormalization of 
the theory on S 4 ^ e . 


2 The Wilson-Fisher fixed points in curved space 

The renormalization of interacting scalar field theory in d = 4 — e on a curved space was studied 
in 38,39) for the single scalar theory, and in [8,40 for the more general multicomponent theory 
(see Appendix [A]). In the case of N real scalar fields with O(N) invariant interaction, the full bare 
action on a curved manifold is 

S = J d d x^g ((<9 m </>o) 2 + ^ d _ 1} Wo)i + + -Ty°i?(^o) 2 + a^W 2 + boE + cqH 2 ^ , 

( 2 . 1 ) 


where 1Z is the Ricci scalar, H = and 


^ + {d _ 2 { d _ x) 

e = n llvXp n l,vXp - + n 2 . 


( 2 . 2 ) 


Here W 2 is the square of the Weyl tensor, and E is the Gauss-Bonnet term which is the topological 
Euler density in d = 4 (it is not topological in d = 4 — e, but we will still refer to it as the Euler 


density in what follows). In (2.1), we have separated an arbitrary coupling r?o to the scalar curvature 
from the conformal coupling term The action (2.1) contains all possible terms which are 

marginal in d = 4, which have to be included in order to consistently renormalize the theory in 
d = 4 — e. Note that in d = 4 — e the coupling Ao has dimension e, a o, bo, Co dimension —e, and 770 is 
dimensionless in any d. We have omitted a mass term for the scalar field since we will be interested 
in the conformal theory, and the mass can be consistently set to zero in dimensional regularization. 

Since the theory is renormalizable, the divergences in the free energy F = — log Z on a general 
manifold can be cancelled by expressing the bare couplings Ao, ao, bo, Co, rjo in terms of renormalized 
ones A, a, b, c, rj. The renormalization of the quartic coupling is well known from flat space and 


reads 


^0 — /T ( A + 


{N + 8) 
8ir 2 e 


A 2 + 


/ (N + 8) 2 
\ 647r 4 e 2 


3(3IV + 14) A 
128vr 4 e ) 



(2.3) 
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Requiring = 0 this yields the beta function 


43 


. . JV + 8 2 3(3IV + 14) 3 (33A ^ 2 + 480iVC(3) + 922iV + 2112C(3) + 2960) 4 

PA = -e\-\ — ^ 0 A- ——- A -A H- R -A + .... 


87T 2 


64-7T 4 


40967T 6 


(2.4) 


The renormalization of the curvature couplings takes the form 38,39 


a 0 = p 6 (a + L a ), 
bo = P e {b + Lb ), 
cq = /i e (c + L c + rjL K + rj~L\) 
Vo = {v + L^Zi 1 , 


(2.5) 


where a,b,c,rj are dimensionless renormalized parameters, and L a , Lb, L c , L K , L\, L , 7 are functions 
of A alone with pure poles in namely 






( 2 . 6 ) 


i=l 


2=1 


and similarly for c,rj,n, A. In the case of the quartic 0(N ) theory (2.1), these functions are given 
explicitly by 

«0i + ™ A2 + o(A 3) = ^ + + ^i)v + 0(A 6 ), Lc = 51V + 0(A 6 ), 

e e e e \ e z e J e 

(2.7) 


where the values of the known coefficients are [ 38—40 


N 


aoi = 


021 = — 


120(4vr ) 2 ’ 

boi = Vo j b,n = — 


N(N + 2) 

216(4tt) 6 ’ 

7V(JV + 2)(31V + 14) 


C51 = 


360(4t r) 2 ’ 240(4vr) 10 

N(N + 2)(N + 8)(3N + U) 


h _ (iV + 8 ) 

452 “ ~[2^ bu ' 


(2.8) 


360(4vr; 
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We will be able to determine the coefficient 651 by a direct perturbative calculation on the sphere 
in the next section, and the result is given in eq. (3.9) The functions L K and L\ will not play an 
important role in what follows, but they have the structure 


La = + N [ N _t 2) ~2 + °( A2 ) > L - = °( A4 ) 


2 ( 47 t ) 2 e ( 47 t ) 4 e 


(2.9) 


6 We use a minimal subtraction scheme. 

'Note that, when working on the sphere, we cannot distinguish the divergences associated with the Euler term 

from those associated with the 77 2 term. Therefore, in order to determine 651, we have to assume the result for C51 

obtained in 139,40]. 
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Finally, in the p renormalization, Z 2 is the Z-factor in the flat space renormalization of the p 2 
operator, which is related to its anomalous dimension by |43 


( 2 . 10 ) 


and the function L v starts at order A 4 


log ^ A - + o(x 3 ) 


8 ir 2 


128?r 4 


38 


39 


L v = mi ^- + 0(\ 6 ). 


( 2 . 11 ) 


The coefficient 7/41 may be read off from 139 40 , but we will not need its explicit form here. Note 


that the A independent terms in L a , L b and L A , which are proportional to N, are those required to 
cancel the divergencies of the free field determinant, Ff ree = ^logdet V 2 + 4 ^_ 2 -^ 7 Z + 
which may be extracted for an arbitrary manifold using heat kernel methods. 

After expressing the bare couplings in terms of renormalized ones, the free energy is finite in the 
limit e —> 0 for any values of the renormalized parameters A, a, b, c, p. The renormalization process 
implies that the curvature parameters a, b, c, 77 acquire a scale dependence, just like A, and have 


their own beta functions. Requiring that the bare couplings are independent of p, from (2.5) one 
obtains the beta functions 

n da * /\\ 

Pa = p— = ea + Pa{A) , 

dp 

Pb = = e6 + /3 b (A), 

dp 

A c 

Pc = p-f- = ec + P C (X) + r]P K (X) + r/ 2 P\(X ), 

dp 

r, dr) n f w 

Pv = = r ll<tP + Pv(X), 


( 2 . 12 ) 


where the hatted beta functions depend on A only (and not on e) and are given by 


Pa — — ^a — At^T^ ~ a oi + 3 a 2 iA" + ... , 

Pb = (e - Px^A L b = (1 + A^) 4 1} = boi + 56 4 iA 4 + 6651 A 5 + ... , 
Pc = (e - Px-^j L c - p,L K = (1 + A ^)L« = 6 c 5 iA 5 + ... , 

Pr, = (jV - P^^A L V = X §x L ? = 4r ?4lA 4 + • • • , 

Pa = (e - Px-^ ~ 2 1(j> A La = (1 + A -^)L^ = 0( A 0 ), 

Pn = (e - Px-^ - 7 ^) L k - 2 p v L A = (1 + A^) 4 X) = ^(A 4 ) , 


(2.13) 


7 











where L^\ L^\ ..., denote the residue of the simple pole as defined in (2.6). As usual, the beta 


functions depend only on the coefficients of the simple pole, but there is a tower of consistency 
conditions (’t Hooft identities) relating higher poles to lower ones. 

Since the bare couplings are independent of the renormalization scale, the free energy satisfies 
the identity 

(2.14) 


o, fro, c 0 , vo) = 0, 

dfj, 

which implies that when F is expressed in terms of the renormalized couplings, it obeys the Callan- 
Symanzik equation 


d d d d d d 
^ + Pvir. + 7 + Pb^_ + Pc 


dfi 


l dX 


dv 


da 


’db 


dc 


F = 0. 


(2.15) 


Note that the dependence of F on the renormalized parameters is simply linear in a, b, c, since the 
corresponding curvature terms are held independent additive terms, therefore we may split 


F = F^X, v, m ) + At e a / dPx^fgW 2 + /i e b d d Xy/gE + /j, e c / d d XyfgH 2 , 


(2.16) 


where F ( p(X,v,ia) is finite for any A, 77 . 

So far we have been working on a general curved manifold. Let us now specialize to the case 
of a round sphere S d of radius R. In this case, we have 'TZp V \ p 'Ry ,vXp = 2d(d — 1 )/ i ? 4 , TZ^TZ^ = 
d(d — 1 ) 2 / R 4 and 7 Z = d(d — 1)/R 4 , therefore 


W 2 = 0, E = d(d-l)(d-2)(d-3)/R 4 , H = d/R 2 


(2.17) 


On the sphere, the radius dependence of the free energy is constrained to be F = F{ A, a, b, c, r/, /j,R), 
and so we see that 

r)F rll? r 

(2.18) 




where we have used the conventions T^ v = Therefore from the Callan-Symanzik equation 


(2.15) we can deduce the relation 


^ d d Xy /g(Tp =P^+Pr^+ F~ e Pb d d XyfgE + J d d x^H 2 . 


(2.19) 


Note that ^ and ^ are related respectively to the integrated one-point functions of the renor¬ 
malized cj) 4 and IZcj) 2 operators. In fact, an unintegrated version of this equation (up to equations 
of motion and total derivative terms, and including the Weyl square term as well) also holds as 
an operator equation on a general manifold. [^] In order to have a scale invariant theory in curved 
space in d = 4 — e, we see therefore that we should set to zero not only the usual beta function for 


’See for instance eq. 10.1 of 39 . 







A, but also all curvature coupling beta function^] 


/ 5 a = / 5 « = /3b = Pc = Pr, = 0 . 


( 2 . 20 ) 


Using the expressions for the beta functions (2.4), (2.12), these conditions can be solved in d = 4 — e 
to give 


A — A*, 


CL* — 


/5a ( A*) 


6 * = — 


/5b(A* 


c* — 


/5 C (A*) + ??*/5 k (A*) + 7*/ 5 a (A* 


V* = 


Pr](K) 

70 2 (A*)’ 


( 2 . 21 ) 


where A* = 


8 tt 2 

iV+8 


e + ... is the value solving /3\ (A*) = 0 which defines the Wilson-Fisher fixed point 


in flat space. More explicitly, using (2.8) and (2.13) 


a* = -3a 2 i^ + 0(e 2 ), 

A 5 

c* 6 r.-,i -* • Ob 5 ). i]* = ~ 


&* = -—- 56 4 i— - 6651^ + 0 (e 5 ) 

e e e 

+ 0 (e 4 } _ 

7</> 2 (A*) 


( 2 . 22 ) 


Tuning all couplings to their fixed point values, we then get the sphere free energy at the IR fixed 
point 

Fir (d) = F(A*, a*, 6 *, c*, r/*, /xR). (2.23) 

We emphasize that the resulting Fir(c/) is just a function of d = 4 — e, the radius and /r dependence 
drop out as a consequence of setting all beta functions to zero. We will verify this explicitly by our 
perturbative calculation below. 


Eq. (2.12) is our definition of the Wilson-Fisher fixed point on a curved space. Note that it is 


important that we work in d = 4 — e with a non-zero e. If we take the e —> 0 limit first, then the 
curvature beta functions become independent of a, b, c and it is not possible to set them all to zero. 
Indeed, if we take d —> 4 first, we get (since A* = 0 in d = 4): 


[ d*Xy/g(Tli) = 6 01 [ d^gE = -^- 
ls 4 Js 4 


(2.24) 


where we used J s4 d^x^fgE = 64vr 2 . This is the expected conformal a-anomaly of N free massless 
scalars 9 10 A related fact is that if we dial to the fixed point (2.12) in d = 4 — e, and then take 


0, the free energy diverges: it has a simple pole whose residue is related to the d = 4 anomaly. 


It follows that the function Tj 74 (d) = — sin(^)lq 74 (d) has a smooth limit as d —> 4 proportional to 


9 Of course, from a calculation of F on S d we are not sensitive to f5 a , but we have included it here since it would 
play a role more generally on a not conformally flat space. 

10 On a general manifold the term proportional to the square of the Weyl tensor aoi f s4 d 4 x^/gW 2 will of course 
also appear on the right hand side. 
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the a-anomaly. Indeed, since the first term in (2.16) is finite as e —> 0, we have 


lirn F ni (d) = -^/3&(A* = 0) f (Px^/gE 

4 2 J s 4 

- < 2 - 25 ) 

While here we are considering the scalar field theory which has a trivial free fixed point in d = 4, 
the fact that lim^-^ Fm(d) reproduces the 4d a-anomaly coefficient is general and should apply 
also to cases where there is a non-trivial fixed point in d = 4. 


3 Perturbative calculation of F for the O(N) model in d = 4 — e 


In this section we perform a calculation of F at the IR fixed point of the 0(N) model in d = 4 — e 
to order e 4 , i.e. we determine F to order e ' 5 (the case of most general quartic coupling is worked out 
in Appendix [A]) . Because of the structure of (2.22), this requires knowing the beta functions for the 
curvature terms to order A 5 . Let us first argue that to this order we can neglect the effect of the 
renormalization of the conformal coupling parameter r]. Note that once F is expressed in terms of 
renormalized parameters, it is finite for any value of 77 . Since the curvature counterterms Lb, L c in 


eq. (2.5) are independent of 77 , we can fix them by working at 77 = 0. This means that we can set the 


bare parameter to ? 7 o = = 7741 ^ + 0 (A 5 ). The leading effect of this parameter is a divergence 

in F of order coming from a two point function diagram ( G 2 in figure [lj) with a 770 insertion in 
one of the propagators (the 0 -point and 1 -point functions contribute to higher orders, respectively 
order F and £). Hence, it does not affect the calculation of Lb and L c to the order A 5 that we 
consider below. After the ? 7 -independent part of the counterterms is fixed, we can reintroduce the 


77 dependence perturbatively. Recall that at the fixed point we set 77 = 77 * = 0(e 3 ), see eq. (2.22). 
The leading contribution to F due to the 77 * is then a finite term of order 77 *A* ~ 0(e 5 ), coming 
from the two point function with a “mass” insertion on one propagator. In addition, the value of 


c* at the fixed point is shifted by the 77 dependent terms in (2.21) only starting at order e 6 . We 


conclude that, to the order we will be working on, we can ignore renormalization of the conformal 
coupling parameter 77 and work with the action 


S — J d d xy/g {(d^b ) 2 + 4^ __ ^^(^o) 2 ) + (^o^o) 2 + ^0 E + cqH 2 

Introducing the integrated 77 -point functions of the free theory 

/ n 

n d ?x i ^ i (4{x 1 )...4(x n ))o im , 

i= 1 


(3.1) 


(3.2) 


10 









where (f >g = (^q^q) 2 , the free energy up to fifth order is given by 


F - F„„ = + F°_G 3 -^04 + + J d d x ^g(b 0 E + off 2 ), (3.3) 

where Ff ree is the contribution of the free field determinant 11 

Ff ree = N log det ( -V 2 + J ~ 2 n) = —N f duusimru F (2 + n ) r (2 ~ u ) = + 0 ( e °). 

V 4(d-l) J J 0 sin(f)r(l + d) 90e 1 7 

(3.4) 

One can obtain the 

e-expansion of Ff ree to any desired order. 

Using the two-point function 


The divergence of this determinant fixes 601 = — 3 60 (^ 7r )2 as given in (2.8 


{4>\x)(P(y))o = 




s(x,y) d 2 

where s(x,y ) is the chordal distance we find 


_ r(f -!) 

* 4vr d / 2 


(3.5) 


G 2 = 8N(N + 2)C 4 4> I 2 (2(d-2)), 

G 3 = 64 N(N + 8 )(N + 2)Cj/ 3 (2d - 4), 

G 4 = 32(4!)1V(1V + 2 )((N 2 + 6N + 20)G'f ) + 8{N + 2 )G% ) + 4(5iV + 22)G^) , (3.6) 

G 5 = 2 12 • 3N(N + 2) (20(N + 2 )(N + 8)g£ ] + 10(3iV 2 + 22N + 56)G^ 5 

+ (IV 3 + 8 N 2 + 24iV + 48)Gf } + 10 (N + 2) 2 Gf ] + 20 (N 2 + 201V + 60) G ( 5 5) + 8(5N + 22)G^ 6) ) , 


where all diagrams are represented in Figure [lj 



Figure 1: Diagrams contributing to the free energy in c/> 4 -theory up to fifth order. Each line 
represents the sphere propogator ((/>(x)0(y)) = C^/s(x, y ) d ~ 2 . Symmetry factors are not included 
in these graphs. 

The integrals / 2 (A) and / 3 (A) are given in Appendix |Bj also using the method described in this 
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appendix we find for the four point functions 


Gf ] = 


3(47r) 8 e 2 


+ 


24(7+log(47r_R 2 ))+73 + 5 (24(7+log(47ri? 2 ))+73) “+360 tt 2 +6713 


9(47t) 8 € 


1080 (47t) 8 


+ ®( e ) 


1 4(7+log(47r_R 2 ))+15 | ^ ^ 

20(47r) 8 e 40(4ir) 8 , 

2 | 120(7+log(47r_R 2 ))+383 | (l20(7+log(47rR 2 ))+383) “+1800 tt 2 +35041 ( /n( 

k 3(4 tt) 8 <e 2 ^ 90(4tt) 8 £ ' 10800(4tt) 8 ^ 

and for the five point functions we get 

5 5(30(7+log(47ri? 2 ))+109) | /nf 

33( 47r )io e 2 2 2 3 4 (47r) 10 e h U V e 1 ’ 

16 . 4(l0(7+log(47rR 2 ))+3l) 3(l0(7+log(47rJ7 2 ))+3l) +30 tt 2 +641 q 

9(47r) 10 e 3 ' 9(47r) 10 e 2 ' 2-3 3 (47r) 10 e ' ' / J 

40 ( 20(l5(7+log(47r_R 2 ))+44) lo(l5(7+log(47ri? 2 ))+44)“+2257r 2 +4612 


cf = < 


9(47r) 10 e 3 3 3 (47r) 10 e 2 2-3 4 (47r) 10 e ) > 

27(47r) 10 5 

8 , 405(7+log(4-7r R 2 ))+49) ( 2(l5(7+log(47ri? 2 ))+49) +45 tt 2 +1014 | ^ 

9 ( 47 r) 10 e 3 4 3 3 ( 47r )lO e 2 I 2-3 4 (4 7 r) 1 °e ^ 1 ’ 

8 C(3) 1 


+ Off °i 
9(47r)lO e + 1 


(3.7) 


(3.8) 


Combining all results in (3.3), and expressing Ao,5o> c o 411 terms of renormalized couplings using 


(2.3) and (2.5)-(2.7)-(2.8), we find that all poles cancel provided we set the unknown 651 coefficient 
in the Euler term to 


N(N + 2) (3(24 + 7N) + 4(57V + 22)C(3)) 
90(4tt) 12 


(3.9) 


This allows to find the Euler beta function /3b to order A 5 , given in eq. (1.7), and hence determine 
the fixed point value 6 * to order e 4 by (2.22). Using the fixed point coupling 


A* 


8vr 2 ^ | 24(31V + 14)tt 2 2 
JV + 8 C+ (N + 8) 3 6 

7 r 2 (33 N 3 - 11(UV 2 + 96(1V + 8)(51V + 22)C(3) - 17601V - 4544) 3 ( 

(N + 8) 5 6 + 


(3.10) 


In this equation NF S (4 — e) = — sin(4^)Ff ree . Note that the dependence on fiR has dropped out, 
consistently with the conformal invariance. As a partial check of this result, we note that at large 
N it yields 


Fir - NF S U - e) =-e 3 -e 4 + 

K v J ' 576 6912 


7T 3 137T 4 


/ 7T 3 6477T A 

V13824 ~~ 414720 ) 


e 5 + 0(e«))+0(i-\, (3.11) 


which precisely agrees with the result derived in [l] using the Hubbard-Stratonovich transformation; 
it follows from (1.5) with A = d — 2. 
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We can also obtain a result for F in d = 4 at a generic value of the renormalized coupling constant 
A. Taking the limit e —> 0 of ( |3.3| ) (after expressing all bare couplings in terms of renormalized 
ones), we find 


Id=4 d=4 _ N(N + 2)A 2 N(N + 2) (N + 8) (5 + 2 ( 7 + \og{^^R 2 ))) A 3 iV(iV + 2)A 4 
free 72(4vr) 4 + 72(4 tt) 6 1440(4vr) 8 

x (l5 (N + 8) 2 (5 + 2 ( 7 + log(4vr^ 2 J R 2 ))) 2 + 90(3IV + 14) (5 + 2 ( 7 + log(4vr^ 2 ii 2 ))) 

+ 631V 2 + 14641V + 6272) + ^tt 2 (36 + 2c). (3.12) 

' O 


For the integrated trace of the energy-momentum tensor due to interactions we get 


d 


- aXyfg(T£) = R q^( F - F hee ) 

_ N(N + 2 ){N + 8 )A 3 N(N + 2)(2 (N + 8) 2 ( 7 + log( 47 r/i 2 i? 2 )) + (5 N + 89)1V + 362)A 4 

“ 18(4vr) 6 12(4vr ) 8 ' 

(3.13) 


As a consistency check of this result, we note that it satisfies the Callan-Symanzik equation (2.15) 
to order A 4 (without the f3 ^ term that we have not included here). 


3.1 Pade approximation of F 


For any quantity f(d) known in the e 
Pade approximant 

Pade[ mn ] (e) 


= 4 — d expansion up to a given order, we can construct a 

_ Ao + A\e + A 2 € 2 + ... + A m e m , . 

1 + Bie + B 2 e 2 + ... + B n e n 


where the coefficients Aj, Bi are hxed by requiring that the expansion at small e agrees with the 
known terms in /(4 — e). If a quantity is known in the e-expansion to order k o, one can construct 
Pade approximants of total order m + n = fco. 

An improved version of the Pade approximant may be obtained if one knows a few terms of the 
e-expansion of the quantity f(d) near two values of dimension: 


f(d\ + e) —do + cqe + a 2 ^ + ..., (3.15) 

f(d 2 - e) = 6 0 + h e + b 2 e 2 + ... . (3.16) 


Using the known expansions near the two boundaries, one can fix more coefficients in (3.14) and 
hence use Pade approximants of higher order, which should provide a better resummation of the 
quantity f(d). 

For the Ising model (N = 1) in 2 < d < 4, we know the expansion of F in d = 4 — e to order 
e 5 , eq. (1.6). In addition, we know that for d = 2 the IR fixed point of the quartic scalar field 
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theory should coincide with the 2d Ising CFT, which is known to have central charge c = 1/2. This 
corresponds to the value = 2) = n/12. In other words, we have: 


-Flsing {d) 



at 

e) T Tint at 


d = 2 , 
d = 4 — e, 


(3.17) 


where 


F s (4-e)+F in t = -^-+0.0205991e+0.0136429e 2 +0.00670643e 3 +0.00264883e 4 +0.000927589e 5 +0(e 6 ). 
180 

(3.18) 

Imposing the d = 2 value as boundary condition then allows us to use Pade[ m n ] approximants with 
total order m + n = 6 . All approximants Pade[ lj5 ], Pade[ 2j4 ], Pade[ 3]3 ], Pade[ 4i2 ]and Pade^q appear 
to give very close results in the whole range 2 < d < 4. In d = 3, we get 


0.977973 

with 

Pade[i j5] 

0.97383 

with 

Pade [ 2 i 4 ] 

0.977106 

with 

Pade [3i3] 

0.977248 

with 

Pade [4)2 ] 

0.975379 

with 

Pade [5)1 ] 


(3.19) 


Notice that all these values are lower than one, consistently with the T-theorem in d = 3. Given 
that these approximants are very close to each other, it seems reasonable to average over them. In 
this way, we obtain the estimate in d = 3: 


« 0.976 . (3.20) 

Fs 


Using instead Pade approximants of total order m + n = 5, which are insensitive to the order e 5 
in F in d = 4 — e, we find a very close estimate Fj s [ ng /F s « 0.974, which indicates that the Pade 
resummation appears to be reliable^] In FigureJ^J we plot Fising normalized by the free scalar result 
F s , comparing the result of the two-sided Pade approximants and the unresummed e-expansion. 

In the general O(N) case for N > 2, the IR fixed point of the quartic theory is known to have an 
equivalent description as the UV fixed point of the non-linear 0(N) sigma model, which is weakly 
coupled in d = 2 + e. For d = 2, one has a free theory of N — 1 massless scalar, which has c = N — 1, 
or F = ( N — 1)7 t/6 . The same d = 2 value applies to the case N = 2, where one has a compact 
scalar with c = 1 (in this case, though, the behavior of the model in d = 2 + e is not known). 


11 Using the unresummed e-expansion (3.18 ) and setting e = 1, one obtains instead Fi s i ng /F s 
and Fi s i r ,g/F 3 S3 0.971 to order e 5 . 


0.957 to order e 4 l], 
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^Ising/^s 



Figure 2: F/F s for Ising model in 2 < d < 4. The solid line is the result obtained by averaging 
Pade approximants Pade[ m!n ] with m + n = 6 and fixed d = 2 boundary condition. The dashed line 
corresponds to the e expansion result (3.18) without resummation. 


For N > 2, the expansion of Fq(N) in d = 2 + e is expected to take the form 

-Fo(at)( 2 + e) = (./V — 1)F S (2 + e) + a\e + a 2 + ... . (3.21) 

where the coefficients 01 , 02 , ... may in principle be fixed by doing a perturbative calculation on 
S 2+e . This suggests that it may be more natural to apply the Pade resummation procedure not to 
the full F 0{N ) = NF S + Fi n t, but to the quantity 


fid) = F 0{n) -(N- 1 )F S = F S + , 


(3.22) 


since this is expected to have an expansion in integer powers of e both near d = 4 and near d = 2 
(at least for N > 2)J^] The exact value at d = 2 is /(2) = 0, and in d = 4 — e we can read off the 


expansion up to order e 5 from (1.6). Once the two-sided Pade approximant of f(d) is obtained, we 


can reconstruct the full by adding the contribution of N — 1 free scalars, using eq. (1.4). 


The d = 3 estimates of Fo(n) f° r a f ew values of N obtained by using the Pade[ 4 2 ] two-sided 
approximant on (3.22) are given in Table [lj For comparison, we also list the values obtained from 
ordinary Pade approximants (without d = 2 constraint) carried out on the full Fo(n) ■ 

As a test of our results, we note that in the large N limit, our Pade extrapolation with d = 2 
boundary condition yields, in d = 3: 


F 0(n) ~ 0.063807N - 0.015398 + 


0.016974 

N 


+ 0(1/N 2 


(3.23) 


2 Note that the expansion of the free scalar answer F a by itself contains log(e) terms in d = 2 + e. 
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The term of order N° is within one percent of the exact value — £(3)/(87r 2 ) « —0.0152242, which 


may be obtained from eq. (1.5) by setting A = 1, d = 3. We have verified a similar good agreement 


with (1.5) over the whole range 2 < d < 4. The term of order 1/N is not known exactly in d = 3, 
and it would be interesting to compute it directly by large N methods and compare with our 


estimate (3.23). 


N 

1 

2 

3 

4 

5 

10 

Pade (b.c.) 

0.9763 

0.9740 

0.9748 

0.9759 

0.9772 

0.9831 

Pade (free) 

0.9779 

0.9765 

0.9762 

0.9766 

0.9773 

0.9820 


Table 1: List of Pade extrapolations for Fq/ n \/(NF s ) in d = 3. The first line corresponds to the 
constrained Pade approximants with d = 2 boundary condition, and the second line to ordinary 
Pade approximants. For N = 1 (b.c.) the entry is the average of several Pade approximants, eq. 
(3.20). For N > 2, the constrained Pade is carried out on the quantity f(d) defined in (3.22). as 


explained in the text, 
using Pade[ 3)2 ]. 


The “free” Pade approximant is carried out on the full Fq(n)i eq. (1-6), 


We note that that the ratio Fq^ n ^/(NF s ) in d = 3, given in Table |l| has a minimum around 
N = 2 using the constrained two-sided Pade, while it has a minimum close to N = 3 using the 
unconstrained Padej^] We observe that a similar behavior, with a mimimum close to N = 2, was 
found in the conformal bootstrap approach 44 for the ratio c'^ 1 °^ N ' > /(Nc^). The estimates for F 


given in Table [T] imply in particular that, in d = 3 


NFi s i n a- > F ( 


0(N) 


for N < 4 . 


(3.24) 


As a consequence of the d = 3 ^-theorem, this suggests that for N > 5 one cannot flow from N 
copies of the 3d Ising CFT to the 0(N) invariant fixed point, while such a flow is possible for 


N < 4. In section 4.4 we will see that a similar behavior can be deduced in d = 4 — e by studying 


an 0 (N ) breaking deformation of the quartic scalar field theory. 

As an additional check of our results, note that by adding a relevant operator (the mass term 
with a negative m?) it is possible to make the CFT flow from the O(N) invariant fixed point to the 
phase where the 0(N ) symmetry is broken to 0(N — 1). In this phase there are N — 1 Goldstone 
bosons, and the sphere free energy far in the IR is 




The F-theorem then requires that 


Fo(n) > F$b 


(3.25) 


(3.26) 


li The latter result is closer to what was observed in 'l^ using the unresummed e expansion. 
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For large N, this inequality was shown to hold in 16 . Using our results obtained from the e 


expansion and Pade extrapolation, we have verified that it in fact holds for all N. 

4 F in Conformal Perturbation Theory 

Let us consider a conformal field theory CFTo in dimension d perturbed by a weakly relevant 
primary scalar operator 

S g = ScfTo + 9b j d d xO(x ). (4.1) 

Here g b is a bare coupling constant, and the bare operator O has dimension A = d — e in the 
unperturbed CFT, with e< 1. The operator O has flat space two- and three-point functions given 
by (in the unperturbed theory) 


(O(x)O(y)) 0 = 


C 2 


\x — y | 2A ’ 
(0(x)0{y)0{z)) o = 


C 3 


(4.2) 


x — y\ A \y — z\ A \z — x\ 


The renormalization of the coupling constant can be deduced from a perturbative calculation of 
the two-point function in flat space. We have 


(0(x)0(y)) g = (O(x)O(y)) 0 -g b d?z(O(x)O(y)O(z)) 0 + G(gl) 


C 2 


\x - y | 


2A 


9b- 


C 3 


x - y\ L 


d d z 


| y — z\ A \z — x\ 


+ 0(g 2 b ). 


(4.3) 


We can use the integral 


d d z 7 


n l r(|-a)r(f-/3)r(a + /3-0 


\x — z\ 2a \y — z\ 2 P \ x _ y\2(a+p-%) r (a) T(P)T(d-a-P) 
so we hnd, after plugging A = d — e, 


(4.4) 


(0(x)0(y)) g = 


In the limit e —> 0, we get 


C 2 


c 3 \x-yf ^r(|) 2 r(|-e) 


x — y| 2 ( d e ) \x — y| 2 ( rf e ) p ( ^ 2 £ ) 2 r (c) 


tt2 ^^ + 0(gl). 


(4.5) 


{0(x)0(y)) g = 


C 2 


x _ y |2(d-e) gb \ x _ y| 2 (d-e) p (|) 


a . . 

C 3 47T2 \x — y| e 


+ 0{gl). 


(4.6) 
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We can introduce a dimensionless renormalized coupling g by 


9b = 9 [ 9 + zi — + • • 


(4.7) 


where g is the renormalization scale. This corresponds to the renormalization of the local operator 
0{x) = Z a O ren (x) , Zo 1 = d{f "^ 9b) =1 + 2 Zl g - + 0(g 2 ) . (4.8) 


—™=l + 2 Zl ?-+ 0 (g 2 ). 
dg e 

We can fix z\ by requiring that the two-point function of the renormalized operator is finite in the 


limit e —> 0. From (4.6) we find 


(Oren(cc)O r en(y))(ji — 


x — y\ 2 ( d ~ e ) 


„ . fl. „ 47T2 (u\x — v\Y 

C 2 (l + 4^i—) - gCz-j-^ -— + 0(g 2 ) 

e F(f) e 


This gives 


zi = 


_ a 

C 3 7 T 2 


d\ ' 


c 2 r© 


From (4.7) we then find the beta function 


(4.9) 


(4.10) 


Pio) — ~ e 9 + jr9~ + o(g 3 ) ■ 

r (i) c 2 


Then there is a IR stable fixed point given by 

r (1) C 2 


g* - 


7 T 2 


d C 3 


e + 0 (e 2 ). 


The dimension of the renormalized operator O ren (x ) at the fixed point is given by 
Air = d — e + /3(g)— log Zo\ g = gt = d + (3'(g*) = d + e + 0(e 2 ). 


(4.11) 


(4.12) 


(4.13) 


Note that the operator has become irrelevant at the IR fixed point. Setting g = g* in (4.9), we find 
that the two point function of the renormalized operator takes the form, to leading order 


(Oren(3')O re n(y))lR — fJ> 


_ 2(d—e) 


(C 2 + 0 (e 2 )) 


(H\x — y |) 2 (d+e+0(e 2 )) 


(4.14) 


Note that the leading correction to the two-point function normalization is expected to be of order 


e 2 . 
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4.1 The sphere free-energy 


We now conformally map the CFT (4.1) to the sphere S d . The two and three point functions of the 

replacemei 

2R\x-y\ 


bare operator are the same as in (4.2) with the replacements \x — y\ —> s(x,y), etc., where s(x,y) 
is the chordal distance 


s(x,y) = 


(4.15) 


(l +x 2) 1 /2(l +y 2)l/2 • 

j.2)2 • 

Working in perturbation theory, we can compute the change in the free energy in the perturbed 


and the metric of the sphere is ds 2 = 


CFT to be 16 


9b, 


9b. 


5F = F-F 0 = ~^C 2 I 2 {d - e) + ^ C 3 I 3 (d - e) + 0{g A b ) 


^9 2 


1 


C 2 I 2 (d - e) + g le g A C 3 -y e I 3 (d - e) - 


d 

7T 2 


(4.16) 


^(|) 


I 2 (d-e) +0(g A ), 


where in the second line we have used the relation (4.7),(4.10) between bare and renormalized 
coupling, and we have used the two and three point integrals on the sphere / 2 (A) and / 3 (A), which 
are given in Appendix |Bj Expanding at small e with fixed d, we find 

W V (— ^) 

(4.17) 


W - e) = 21 + 0(e 2 ), J 3 (d - e) = — ^ ^ + <D{e) 


r m 


T(d) 


Note that 


h (d — e) 87T2 


h(d-e) eT (|) 
To leading order in e and g we then find, using (|4.17|) 


+ O(e 0 ). 


(4.18) 


21-d^d+^Y (_<A 

r (¥ ) ^ 


1 


2 1 7T 2 C3 3 


2 e9 + 3r (0 c 2 


9 


(4.19) 


Finally, in terms of F = — sin(^-)E, after some simplification of the d-dependent prefactor we 
obtain 


5F = 


2C 2 n 1+d 

r(i + d) 


1 2 . 1 71-2 C 3 3 

~ 2 cg + sf(|yo 9 


Note that to this order we find 


85F 2C 2 ir 1+d 


Pig) 


(4.20) 


(4.21) 


8 g F (1 + d) 

and so F is stationary at the fixed point. To evaluate the change in F from UV to IR, we can plug 


19 




























in the value of the fixed point coupling (4.12), and we get 


' ~ ~ (|) 2 cl , 

SF = FiR _ FuV = ^ 

In odd d this result agrees with |l 6 | , and in even d with the change in a-anomaly computed in |5j. 
For a unitary CFT, C 2 > 0 and C 3 is real. So we find that for all d, F\jy > Fir to leading order 
in conformal perturbation theory. This gives a perturbative proof of the Generalized F-theorem in 
the framework of conformal perturbation theory. 

we may rewrite the leading order change in F as 

2ir 1+d r 9 * 

,F ~W^L MM- < 423 > 

It is natural to ask whether such a relation continues to hold at higher orders, provided we replace 
C 2 by a coupling dependent two-point function coefficient C 2 (<?). 


Note that using (4.21 


4.2 Several coupling constants 

More generally, we can consider a perturbation by several primary operators 

Sg = FcFTo + E ^ ^ ^ ’ ( 4 ‘ 24 ) 

i 

where g l b are bare coupling constants. For simplicity, we take the bare operators O* to have the 
same dimension 

A = d-e, (4.25) 

where as before we assume e< 1 . 

The two and three point functions in f 

(Oi{x)Oj(y)) 0 = 
i(x)Oj(y)O k (z 

Here Co,ij is a symmetric matrix (positive in unitary theories). We could choose a basis of operators 
so that C 2 ) ij = 5ij, but we will not do this here. 


iat space take the form 

^2 ,ij _ 

>A ’ 

Cjjk 


\x — y | 2A ’ 
))o = 


(4.26) 


x — y\ A \y — z\ A \z — x\ 
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Let us consider the calculation of the two-point function in the perturbed CFT. We have 


{0^)0,(y)) g = (OiWOjiyVo-^gb / d d z(O i (x)O j (y)O k (z)) 0 + 0(g%) 


C‘2.i 


\x - y | 


2A 


E 


„k '■'ijk 

9b | | / 

\x - y\ l 


crz 


| y — z\ A \z — x\ 


+ 0(g 2 b ). (4.27) 


Using the integral (4.4), and expanding at small e, we find 

C-2,ij 


(Oi{x)Oj(y))g = 




x - y| 2 ( d " e ) ^ b \x- y| 2 U- £ ) r (|) 


We can introduce the dimensionless renormalized couplings g t by 


gi = in 9 i +Y,zj 


g j g k 


jk 


+ 


jk 


The corresponding local operators renormalization is 


(4.28) 


(4.29) 


Oi(x) = (Z 0 )ijO™(x), 

Ojr(x) = {S i k + 2z* k gi)O i (x). 


(4.30) 


Plugging this into (4.28) and requiring that the renormalized two-point function has no poles as 


e —> 0, we get the condition 


2 C 2 ,ikzLg m + 2 C 2tjk zf m g m - ^c ijm g m = 0 . 


k jm 


d 

47T2 


r (i) 


This is solved by 


4 = c$c ljk 


d 

7T 2 


d\ ’ 


r (S) 


where Oj is the inverse of the matrix C 2 u- The beta functions are then given by 


(4.31) 


(4.32) 


P i = g d ^- = -eg 1 + -^C 2 l C ljk g^g k + 0(g 3 ). 
ay- i ( 2 ) 


d 

7T 2 


(4.33) 


The computation of the sphere free energy proceeds as in the single coupling case described 
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above. In terms of the integrals I 2 and I 3 we have 


9b9b, 


9b9b9b. 


5F = F-F 0 = C 2 tlj I 2 (d - e) + C ijk I 3 (d - e) + 0{g A b ) 

d 

-Wg j C 2 ,ij + l --^C ijk g i gig k + 0 (g A ) 

z 6L UJ 


2 1 -V+ir(-|) 

rm 


(4.34) 


where in the second line we have used the small e expansion of I 2 and I 3 , and we have plugged in 
the relation between bare and renormalized couplings. In terms of F, we have 


5F = 


2 ir 1+d 


r(i + d) 

Analogously to the one-coupling case, we then find 

86 F 2vr 1+d 


d 

+ l^ja\C ijk g l g J g k 


3r (|) 


A/ • 


(4.35) 


(4.36) 


dg l T(l + d) 

To write the final result more explicitly, we note that the vanishing of the beta functions implies 


d 

7 T 2 


0 = ~eglgiC 2 ,ij + =jjrC ijk glglg k 

1 UJ 


and so we find 


5F = Fir — Fuv = —; 


7 r 


l+d 


-C 2 , ijg lgie + 0(e 3 ). 


(4.37) 


(4.38) 


3 r (1 + d) 

In a unitary theory, the matrix C 2 .ij is positive definite and the coupling constants are real, so we 
find as expected Fuv > Fir- 

Note that (4.36) implies that the leading order change in F from UV to IR can be written as 

27T 1+d /F* 


the line integral 


5F = Fir - Fuv = 


Coijpdg 1 , 


(4.39) 


F (1 + d) J 0 ~ ZA3h 

where the integral is along a path connecting the origin in coupling space to the point g * = 


■ ■ •)■ The integral is independent on the choice of path due to (4.36), which implies 

C ^Qgk ~ C ^Qgi • 

4.3 The 0(N ) model in 4 — e dimensions 

Let us consider the scalar field theory with O(N) symmetric interaction in d = 4 — e 


(4.40) 


s = j d d xQ(d M ^) 2 + |(^) 2 ) 


(4.41) 
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We can view the quartic interaction as a perturbation of the form (4.1), where the unperturbed CFT 
is the free theory, and the perturbing operator 0(x) = \(4>i4>i) 2 has dimension A = 2d — 4 = d — e 
in the UV. Then we can formally apply the results of the previous section to the present case. Note 
that in the above calculations we have assumed a CFT defined in dimension d with an operator of 
dimension d — e, and e in principle unrelated to d. Therefore, we are expanding in e first with d 
fixed, and at the end set d = 4 — e (this corresponds to a different order of limits compared to the 
approach in section [3]). In the free theory at dimension d, we have 


(O(x)O(y)), - \n(N + 2) )‘ , 


(4.42) 


(0(x)0(y)0(z)) o = N(N + 8)(JV + 2) 


T(d/2 - 1) 

47T rf / 2 


1 


{\x - y\\y - z\\z - x\ 


\d—e 


From which we can read off, in the notation of the previous section 

4 


C, = !*<* + 2 ) 


C 3 = 


JVT(iV + 8)(JV + 2)( r( ^ /2 1) ) , (4.43) 


so that the conformal perturbation theory /3-function reads 


f3 = -eg + 


( N + 8) r (2 — 1) 2 2 , 


8vrir(f) 


V + W). 


(4.44) 


Note that setting d = 4, this agrees with the well-known result for the beta function in the quartic 
scalar field theory computed in the standard field theory approach, (3 = —eg + ^4^ g 2 + 0 (g 3 ). 

To obtain the leading order change in F on the sphere, we just need to apply the general result 


(4.22). Using (4.43), and setting d = 4 — e, we find 


SF = - 


n N(N + 2) eS + 0(£4) 


576 (IV+ 8) 2 


(4.45) 


This agrees precisely with the leading order term in (1.6), first obtained in [I]. 

In the conformal perturbation theory approach, it is not easy to deduce the next order correction 


in the beta function (4.11), since it depends on the 4-point function of the perturbing operator. 


However, in the present case we know that for d = 4 it should reduce to the known beta function 


computed in the minimal subtraction scheme, eq. (2.4), namely we should have 


0 = ~eg + 


(IV + 8)r (f - 1) 2 /3(3IV + 14) 3 , 


s^r(f) 


\ 64-7T 4 


TO(e) )g A + Otf ). 


(4.46) 


Assuming that the relation of the form (4.23) holds to higher orders, and using the fact that 
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C 2 (g) = C -2 + 0(g 2 ) (see (4.14)), we then find 


SF = 


2 n 1+d 


r (i + d) J 0 


c 2 


(JV + 8)r(^-l) 2 /3(31V+14) 


8-fr (I) 


647T 4 


~ e 9 + 2 - . -- 9‘ - y—^7ZI - I 9 

IT N(N + 2) 3 7 r 1V(1V + 2) (13IV 2 + 3701V + 1588) 4 


576 (IV + 8) 2 


e — 


6912 


(IV + 8) 4 


3 d 5 + 0(e 5 ) 
e 4 + 0 (e 5 ) • 


(4.47) 


Remarkably, this agrees with the result (1.6) obtained by the direct method based on the renor¬ 


malization of the field theory on curved space, provided one includes the effect of the Euler density 
counterterm which affects the order e 4 in F. 


4.4 The 0(N ) model with anisotropic perturbation 

Let us consider the O(N) scalar field theory in d, = 4 — e perturbed by an operator that breaks 
0(N ) to S/v x , which are the symmetries of the IV-dimensional hypercube. This is why it is 


sometimes called the model with “cubic anisotropy” 43 


S= l <iM i(a„« 2 + h( 




(4.48) 


Both operators in the action, (</>*0*) 2 and (]>T </> 4 ), have dimension A = 2d — 4 = d — e in the free 


UV CFT, and so we can apply the results of section 4.2 to this model. In this basis of operators 
the two-point function matrix C 2 ,r/ is not diagonal. While it is not necessary to do so, we find it 
convenient to change the basis to 


S = j d d x + y 


(4.49) 


In terms of the operators 0\{x) = <t>i) and O 2 (x) = \ J2i<j the two-point functions i 

the free theory read 


m 


{Oi(x)Oj(y)) = 


(C 2 ) 


v 


X _ y|2(d e) ’ 


3 /T(d/2 — 1) 

(C2,U - 2 N ( 4 7T d / 2 


IV(JV-l) ( T(4/2-\) 
[ 2)22 2 I 4vr rf / 2 


(4.50) 


(C 2 )l 2 = 0 . 
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The one-loop beta functions computed with the standard dimensional regularization and minimal 


subtraction scheme are known to be [431 (see A.l for the higher order terms) 


/ 3 i = ~egi + 
/? 2 = -eg2 + 


9 g\ + (N- 1 )g\ 

87r 2 

(N + 2 )gl + 6 giff 2 

87r 2 


(4.51) 


Besides the free UV fixed point, there are three IR fixed points, given to the one-loop order by 


87H 


IV + 8 


9i ~ 5-2 — 

5i = ^-e, 5-2=0, 

87T 2 (jV-l) 

Ji 9IV 


* _ 87T 2 

32 3iv e : 


O(N) theory 
IV decoupled Ising 
Anisotropic fixed point 


(4.52) 


The first one is the usual 0(N ) invariant Wilson-Fisher fixed point, the second is the decoupled 
product of N Ising fixed points, and the last one is a new fixed point where the 0(N) symmetry is 
broken by the anisotropic term. The eigenvalues \± of the stability matrix Mjj = ^ 7 , which are 
related to the dimensions of the nearly marginal operators at the fixed point by A± = d + A±, are 
given by 

(4 - N)e 


O(N) : A+ = e, A_ = 
N — Ising : A+ = e, 
Anisotropic : A + = e , 


IV + 8 


A - = -3’ 


(4.53) 


A_ = 


(.N - 4)e 
3 N 


Thus we see that the 0(N) symmetric fixed point is IR stable for N < N c = 4, while the anisotropic 
one is IR stable for N > N c . The IV-Ising fixed point always has one relevant operator which triggers 
a flow to the 0(N ) symmetric fixed point for N < N c or to the anisotropic one for N > N c . Plots 
of the RG flow trajectories for N = 2,4, 6 , 8 are given in Fig. [3} The value N c = 4 is the leading 
order one-loop result. Higher order corrections make N c decrease as we move away from d = 4. In 
d = 3, it is expected that N c ~ 3 143]. 

We now compute the leading order correction to F at each fixed point and check that the 
structure of the RG flows described above is consistent with the Generalized T-theorem. The 


beta functions (4.51) differ from the conformal perturbation theory beta functions (4.33) by the 


d-dependent factors in the quadratic terms. However, this difference does not affect the result for 


F at leading order in e, so we can apply directly our final result (4.38). Using (4.50) and setting 
d = 4 — e, this yields 


Air = Afree — 


jV(3(ffi *) 2 + (iV-l )(g 2 * 
576(47 t) 3 


* A 2 A 


+ 0 (e 4 ) • 


(4.54) 
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It is straightforward to check that this agrees with the result obtained by the approach of 1 , where 
the leading contribution to 5F just comes from the integrated 2-point function on the sphere. 


Using (4.52), we then get for the three fixed points 


F —F N(N + 2)ir 3 p, 4\ 
F °W ~ Ffree 576(A^ + 8) 2 + ° ( } ’ 


Fn-I sing — N. i^lsing — -^free N 
F anisotropic = F Free 


7r 


15552' 3 + 0(£4) ' 

(N + 2)(N - 1)tt 3 


(4.55) 


15552 N 


e 6 + 0(e 4 ). 


Of course, the leading correction to F is negative in all cases, consistently with the fact that we 
can always flow from the UV free theory to any of the fixed points. It is interesting to examine 
the constraints imposed by the Generalized F-theorem on the possible flows between them. For 
instance, from the above result we see that 


F, 


0(N ) ^ F a nisotropic for IV > 4 , 


(4.56) 


which is precisely consistent with the fact that for N > 4 the 0(N ) fixed point has one relevant 
operator which triggers a flow to the anisotropic fixed point, and vice-versa for N < 4. We also 
find that 


IVFi s i ng > F ( 


O(N) 


for N < 10. 


(4.57) 


which forbids flows from the tensor product of N Ising models to the O(N) fixed point for N > 10. 
This is consistent with the structure of the one-loop RG flows depicted in Fig. [3j even though 


the constraint imposed by F appears to be slightly weaker. The result (4.57) just follows from 
the leading order term in F; as discussed in section 


3.1 


including higher orders and doing a Pade 


extrapolation yields (3.24) in d = 3. We also see that 


IV F[ s i n g > Fanisotropic for IV > 2 , 


(4.58) 


which is consistent with the flows between IV-Ising and anisotropic fixed point for N > 4 (again, 
the bound from F appears to be slightly weaker). We also note that at large N the anisotropic 
fixed point gets very close to the IV-Ising point. Indeed we have 


7Tt 1 

IVFj s i n g — Finisol.ropic T T G ( ) ■ 


(4.59) 


The flow from the N decoupled Ising models to the anisotropic fixed point is caused by the “double¬ 
trace” operator Yli<j this is why the correction to F is of order IV 0 . 


As done in eq. (4.47) for the O(N) model, it is interesting to compute the next order correction 
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Figure 3: Structure of the one-loop RG flow trajectories for N = 2,4, 6, 8. The black dot indicates 
the free UV fixed point, the blue one the 0(N) fixed point, the green one the N decoupled Ising 
models, and the red one the anisotropic fixed point. For N < 4, the 0(N ) fixed point is IR stable, 
while for N > 4 the cubic point is. For N = 4, the O(N) and cubic fixed point coincide. 


to F at the anisotropic fixed point using the conformal perturbation theory approach, under the 
assumption that the line integral formula (4.39) holds at higher orders with C- 2 ,ij{g ) = + 0(g 2 ). 


The /3 functions computed in conformal perturbation theory take the form 

Q ( \ | F ( ~ 2 ~) (c\ 2 i ( -\t i 'i 2 \ 51 g\ + 5 (N — 1)5i52 + 4(1V — 1 )g% n ^ 

52 ) = 1 + A 7 Td / 2 {d _ 2) ( 9 ^1 +( N ~ W - ^4 - + °( £ )) > 

£2(51102) - -eg 2 + 47rd/2 ^ _ ^ {( N + %2 + 65152) ^4 (1 + O(e)) . 

(4.60) 

Choosing a piecewise straight path from the origin to the fixed point gl,g 2 , we obtain 
2ir 1+d f 9 2 9 TT 1+d f g i 

Fir - Fuv = — , A / ^2)22^2(0,52)^2 + F7TTFF / (^2)11^1 (51,52)^1 j ( 4 - 61 ) 


r (1 + d) Jo 


F (1 + d) J o 


where the value of 54,52 t° the needed order can be obtained from the /3-functions (4.60). Carrying 
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out the integrations, we find 


Fjr - Fry = - 


ir(N + 2)(JV - 1) 3 n(N - 1)(73 N 3 + 36A 2 + 432JV - 424) 4 


15552A 


-e — 


559872A 3 


e 4 + 0(e 5 ). (4.62) 


This exactly agrees with the result (A.21) obtained using the direct dimensional renormalization 


in curved space, providing a non-trivial consistency check of our approach. 


4.5 The Gross-Neveu model in 2 + e dimensions 

The techniques of conformal perturbation theory may be also applied when the perturbing operator 
is weakly irrelevant and one finds a nearby UV fixed point. As an example, we consider the Gross- 


Neveu model 45 in d = 2 + e: 


S = I d d x (i+ |(^V0 2 ) » 


(4.63) 


where ifj 1 is a collection of N Dirac fermions (we will denote N = JVtrl). The perturbing operator 
0(x ) = has dimension A = 2 (d— 1) in the free IR theory. In d = 2 + e, we have A = d + e, 

and the operator is weakly irrelevant. 

The two- and three-point functions of the perturbing operator 0(x ) in the free theory are 


<0 W 0 to ))^(V-l) hj 


r(|) 


6 


1 


(0(x)0(y)0(z)) = —N(N — 1)(N — 2) . . 

V 27T2 / {\x - y\\y - z\\z - x|) 


(4.64) 


We can formally apply the result (|4.11|) with e —> — e and obtain the beta function 

' N — 2 


o (iv-2)r(§) 2 

P = eg -3- 9 + 


d 

2n 2 


47r 2 


+ G(e) U d + 0( 5 4 ), 


(4.65) 


where we have included the subleading term by requiring that it matches the known beta function 
of the Gross-Neveu model in d = 2 1461. 

Applying (4.23), and again assuming that C- 2 (g) = C 2 + 0(g 2 ), we obtain 


2n 1+d [ 9 * n . . , tvN(N — 1) 3 ttN(N-1)(N-3) 4 /ri . 

5F = ^-- / C 2 p(g)dg = Ao .\ T > e 3 - + 0(e 5 ). 


r (1 + d) Jo 


48(A- 2) 2 


32(A - 2) 3 


(4.66) 


As shown in 47 , this precisely agrees with the direct field theory calculation in curved space in 


d = 2 + e, where the Euler curvature beta function is set to its zero, following the same logic as 
outlined in Section [2j 
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A General quartic theory on the sphere 

In this Appendix we consider, following [8,40 , a general theory with the action 

l/,„ d — 2 l.n U ,1 1 


S = J d d x ^-((a M 2 + + 2 ^ 0 + aoIU 2 + b 0 E + c 0 H 

(A.l) 


where A ^- kl is a symmetric tensor. The /3-functions are given by 


40 


R x -- e \ . 6 a( 2 ) , ° tfl(3,l) iq 0 CT\ 

Pijkl tA ijkl -r (47r)2 Py« f (4^)4 wij'feZ lz, Pijkl ) 


(3,1) 


?(3,2K 


+ 
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2(4vr) ( 


- Ww + + ««*>/$? - /5&?) + -, <A.2) 


where 


Pijkl — ^abij^abkl H - • • • i Pijkl — ^abcd^aijbribed H” • • • j Pijkl — ^abij^ac dk^bcdl + • • • 

Pijkl — ^abef^cdef^acij^bdkl T • • • , Pijkl — ^aede^afij ^bcde^bfkl T • • • , Pijkl — ^abij^acde^befl^cdfk T 
Pijkl ^abij^acdk^befl^cdef T • • ■ j Pijkl ^acdi^aefk^bcdj^befl T • • • , Pijkl ^abciPadf j^dhfk^-jeel T 
Pijkl ^abcd^aikl^jbef ^edef T • • • 


(A.3) 


and the ellipses mean all inequivalent permutations of i,j , k, Z[^] One can recover the usual O(N) 
theory using the tensor 


A ijki 3 A(<%4/ T SikSji T Sudjk). 


(A.4) 


Using /3-functions (A.2) one can find how bare couplings are connected with renormalized ones: 


A ijki ~ P e ( A ijkl + 


(47r) 2 e 


+ -. , 


(A.5) 


For instance p\P kl — ^abij^abkl ^abik^abjl ^abil^abjk 3-11(1 — ^abcd^aijk^lbcd ^abcd^aijl^kbcd H" 


~'ijkl 

^abcd^aikl^jbcd "h ^abcd^ajkl^ibcd- 
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due to cumbersome expressions we only write the first term. For the connected correlators we find 

G 2 = 4!*§C$J 2 (2(d-2)), 

G 3 = 3(4\) 2 t° 3 Cp 3 (2d-A), 

Gi = 12(4!) 2 (9t2 1 Gi 1) + 8tl 2 Gf ] + 36^^) , 

^5 = (4!) 4 (60t5 1 Gg 1 ' ) + 90t® 2 GP + 9*5 3 G® + lO*;^^ + ISOtggGg 5 ' 1 + 2At® & GP) , (A. 6 ) 

where the tensor structures are easily read off from the diagrams in figure [lj 


^2 — ^ijkl^ijkh ^3 — ^ijkl ^ijmn A ki mni 

^41 = ^ijkl^ijmn^mnpr^prkli ^42 = ^ijkl^lmnp^mnpr^ijkn ^43 = ^ijkl^kpmn^mnlr^ijpn 

tbl = Kjkl^kl mn^ijmp^prst^nrsti *52 — A r j kjA /,■ / n i n Apr s n ^prn > t A %j s t j *53 — ^ijkl^klmn^mnpr^prst^ijsti 

*54 — ^ijkl^ijkm^lmnp^prst^nrsti *55 — ^ijkl^lmnp^pnkr^jrst^imsti *56 — ^ijkl^lmnp^pksr^jrnt^imst- 

(A.7) 


Using (A.3) and (A.5) we find 


+0 _ 2 e/. 36*3 / 12(*4 2 — 18 * 43 ) 324(*4i + 2 * 43 ) 

*2 — M I *2 + 71—71 T ( -71—11-T 


(47r) 2 e V (4-7r) 4 e 


(47i) 4 e 2 


18 ( 11*51 + 12 * 5 2 — 168*55 — 96 * 56 C( 3 )) 
(47r) 6 e 


720 (* 5 i — 6*5 2 — 122*55) 2592 (4*5 2 + *53 + 4*55)1 

1 / a \ a o 1 / a \ a q J ' 


(47r) 6 e 2 


(47r) 6 e 3 


0 3 e i 18(*4i + 2 * 43 ) /18(*5i — 6 * 5 2 — 12 * 55 ) 216(4* 52 + *53 + 4 * 55 )' 

*3 ( *3 + -77^:- + ( -- + - / ^4 2 - ) + 


(47r) 2 e 


(47r) 4 e 


(47r) 4 e : 


4 =44 u, U 4(2 /f + * 5s) + ...), 4 = + + ...), 4 = 4'(*« + 24(t5 ?l7 t55) + .. 


(47r) 2 e 


(47r) 2 e 


(47r) 2 e 

(A. 8 ) 


Demanding that the expression for the sphere free energy 


g 2 g 3 

F-F {iee = -^P r 7 + 3 


2 ! • 4 2 3 ! • 4 3 4 ! • 4 4 5 ! • 4 5 


° 4 + + [ d d x^(b 0 E + c 0 H 2 ) (A.9) 


have no poles in e cannot fix the counter terms completely, because E and H 2 are proportional to 
each other on a round sphere. However, using also the result for C 51 from | 8 |, we find 


, - e . &01 . *42 - 18*43 . ( 

6 “ = '‘ l 6+ T + 80(47r) 10 e + 1 ~ 

_ e / 3(*51 — 6*5 2 — 12*55) . 

Co = M c-7 T1 -H 


11*51 + 51 * 5 2 — 186*55 ~ 144 * 56 C( 3 ) 3(*51 — 6*5 2 — 12*55) \ 


40 ( 47 r) 12 e 


+ 


20 ( 4 t r) 12 e 2 


V ")’ 


40 ( 47 r) 12 e 


(A.10) 
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where 6qi = — 360 (^ 7r )2 is the free field contribution determined by (3.4), which just depends on the 


total number of fields. For the curvature /3-functions we have therefore 

^42 — 18*43 3(11*51 + 51*52 — 186*55 “ 144*56^(3)) 


Pb — tb + 601 + 

P c = ec- 


16(47r) 10 
9(*5i — 6*52 — 12 * 55 ) 


20(4?r 


112 


20(4tt) 


12 


(A.11) 


All these results are in agreement with the results for the 0(N ) model and the first terms in 
and p c agree with 1 8 ,40 . The term of order A 5 in pb is, as far as we know, a new result. 


A.l The 0(N ) model with anisotropy 

As an interesting special case, let us consider the 0(N ) model with anisotropy, whose action is 


given in (4.49). In this case our Xijki tensor is given by 


g 2 

Xijki (5 1 9p)$ijkl 4 ._) T T 8il&jk) , 

where dijki = 1 if i = j = k = l and Sijki = 0 in another case. 

The relation between bare and renormalized couplings reads 

£ ( , 9*7i + (3V — l)g| (%\gl + lb(N-l)gigl + (N 2 + N-2)gl 

313 ( 9l + -- + (-64- 

51<7i + 5 (N - 1 )gigl + 4 {N - 1 )sf' 


(A.12) 


128vr 4 e 


+ 


+ 


+ 


^-3 (729 gf + 186(iV - 1 )g\g\ + 16 (N - 1 )(N + 2 )g l9 3 2 + (N - 1 )(N 2 + 9 N- 1 ) ff2 4 )) 
( - 3213,7? - 498(iV - 1 )gjg 2 2 - 10(N - 1 )(N + 38 ) 51 g 3 2 — (N - 1)(53A - 17 )gl 

122887 r 6 e ( 27(145 + 96C(3))ff i + 162(iV " l)g ^ ~ 2(7V “ X )( 13iV " 430 - 192C(3))*7i5 2 3 

+ (N- l)(59iV - 67 + 96C(3))<7?) + 


5127r 6 e : 

1 

3072vr 6 e 2 

1 


c , , 65152 + {N + 2)^2 745*7152 + 9(iV + 2)#ig| + (A 2 + 7iV + 1 ) 5 ! 

52.0 =5 52 H--t- 


87r 2 e 

155i52 + 365152 + 9(iV - 1 ) 5 ! 


647r 4 e 2 


1287 r 4 e 


+ 


+ 


+ 


5127r 6 e : 

1 

1024vr 6 e 2 

1 


(3605^52 + 78 (N + 2)gfg% + 12N(N + 8)515! + ( N + 2)(iV 2 + lOiV — 2 )g, 

( - 384 5 ? 5 2 - (25A + 446)5?5l - 24(8 N + l)g l9 l - 7 (N - l)(3iV + 10)5?) 
122887r 6 e ( 1512 ^ 52 + 9(iV + 138 + 192C(3)5?5 2 2 ) + 24(19iV + 37 + 48C(3))5 i5 2 3 

+ (331V 2 + 457 N - 682 - 768C(3) + 4801YC(3))5?) + ... ). (A.13) 
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and for the /3-functions we find 43 


Pi 


P2 


9ffi + (N - 1 )gl _ 51 gf + 5(1V - l) 9 iff? + 4(1V - 1 )gf l 

87 r 2 647T 4 40967T 6 

+ 162(1V - l)glgl - 2(N - 1)(13JV - 430 - 192C(3)) 5lff | + (N - 1)(59A* 


(27(145 + 96C(3))«rf 
- 67 + 96C(3)) 52 4 ) , 


69152 + (N + 2)g\ 15gfg 2 + 36gig% + 9(N - l)gf 
- £52 H--;n— 7 -H 


9 2 


87r 2 647t 4 40967T 6 

+ 15129 ? + 24(19 N + 37 + 48C(3))9i9| + (33A * 2 + 4571V - 682 + 96(5JV - 8)C(3 ))gf 


+ 138 + 192C(3))9?92 


(A.14) 


The tensor structures defined in the previous section read explicitly 
*2= l -N{?,gi + (N-l)gl), 

*3 = ^N(27g\ + 9 (N - 1 ) 919 ? + (N - 1 )(N + 2)gf), 

*41 = ^(81 gf + 54(IV - \)g\g 2 + 36(1V — 2) (IV — 1)9i9 2 3 + (N - 1 )(N 2 -3 N+ 11 ) 5 ?), 

*42 = ^2 (9ffi + 6(!V — 1 )gfg 2 + {N — 1) 2 9?), 

*43 = ^4 ( 8 I 91 + 18 (N — 1 ) 9 i 9 2 + 24(1V — 1 ) 919 ? + 5 (N — 1) 2 9?) (A.15) 


*51 (39? + (N- 1)9?)(279? + 9 (N - 1)919? + (N 2 + N - 2)9?) , 

*52 (2439? + 108(iV - 1)9?9? + 9(JV - 1)(N + 4 )g 2 g 3 2 + 3(9 N - 7)(N - 1)919? 

+ (N — 1)(31V 2 — 51V + 8)9?) , 

*53 (2439? + 270(1V - 1)9?9? + 90(IV - 2)(IV - l^ 3 + 15(1V - 1)(A^ 2 — 31V + 3)919? 

+ (IV - 1)(!V 3 - 41V 2 + 61V + 12)9?), 

1.54 =p(39i + (IV - 1)92) (39? + (IV - 1)9?) 2 , 

*55 (2439? + 54(IV - 1)9?9? + 72 (N - 1 )g 2 gf + 3 (N - 1)(5!V + 3)919? + (N - 1)(N 2 + 8IV - 12)9?) , 

*56 (8I9 ? + 30(1V - 1)9i9 2 + 15(IV - 1)919? + (N - 1)(5 N - 8)9?). (A.16) 


Writing (note that the coupling constants are absorbed in the coefficients below) 


, -e 7 , 1 Ol31 1 h 41 | / »52 . 051 \ . 

= ^ u +T + T + (_ + T j + ... 


-e / . C 51 

c 0 = M ( c + — 


(A.17) 
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we find, using eq. (]A. 10): 


N (51 gf + 10(IV - 1 )glg 2 2 + 16(IV - 1 )g l9 3 2 + 3(1V - 1 ) 2 g|) 

41 240(4vr) 10 

651 = 540(47r)i2 ( 54 gi( 36 C(3) + 31) + 153(A - 1 )g\gl - (A - 1)(31A - 631 - 720C(3)) ff ^3 
+ (A - 1)(62A + 360C(3) + 169) g^A + (A - 1)(24(5A - 8)C(3) + 95A - 143 )g^j , 
c 5i = 3 gQ(' 4 7 r ^i 2 + 126(A — l)g\g 2 + (A — 1)(5A + ll%)g\g 3 + (A — 1)(35A + l)gig 2 

+ (A- 1) 2 (3A + 10)g^j . (A.18) 


The curvature beta functions are then given by 


/3b = eb + b 0 i + 5641 + 6651 + ... , /3 C = ec + 6 C 51 + ... . (A.19) 


From these we can obtain the values 6* and c* at the IR fixed point /3b = /3 C = /?i = P 2 = 0, as a 
function of the couplings g\, g 2 . At the anisotropic fixed point, the latter are given by 


^ _8^ 2 (iV —l) e+ 8(106- 159A r + 36A 2 + 17A^ 3 )7r 2 e 2 + 


7r 


(a (418304 


9A “ ' 243A" 3 6561A 5 

+ 5848A 3 + 709A 4 - 2592(A - 1)(A 3 + 3A 2 + A - 14)C(3)) - 179776)e 3 + 0(e 4 ) 

8 vr 2 8tt 2 (106 - 125A +19A 2 ) 


- 299734A + 54649A 2 


92 3A 6 


6 2 - 


7r 


■ (a (360640 - 229974A + 41971 A 2 


81A 3 2187A 5 

+ 1955A 3 - 2592(2A 3 - 6 A 2 - 7A + 14)C(3))) - 179776)e 3 + 0(e 4 ). 


(A.20) 


Specializing (A.6) to the present case and using (A.9), we find that the free energy at the anisotropic 


fixed point is: 


Acubic AFf ree 
+ 


(A + 2)(A - l)vr e 3 (A - 1)(73A 3 + 36A 2 + 432A - 424)ir 4 


15552A “ 559872A 3 

--- F ( - 53621A 6 + 14259A 5 + 224610A 4 - 2161280A 3 + 5542944A 2 - 5364672A + 1797760 

302330880A 5 V 

+ 810tt 2 A 4 (A 2 + A - 2) + 10368A(A - 1)(A 4 - 4A 3 + 16A 2 + 6A - 28)C(3)) e 5 + ... . 

(A.21) 


We notice that at 


A r — A rr it — 4 — 2e — ( — — 


5C(3) ^ 


e 2 + 0 (e 3 ) 


(A.22) 


we find F 0 ( N ) = A cu bi c , where Fq is given in (1.6). This is as expected since for this value of A 


the couplings at the O(A) and anisotropic point coincide (compare (3.10) and (A.20)). 
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B Feynman Integrals on the Sphere 


In this appendix we describe a method for computing Feynman integrals on a sphere. We start 
from the well known integrals, which can be calculated explicitly |3,16,48 : 


h{ A) = 


d' i xd d y s fg^ y /(fy 
s(x,y) 2A 

jd™ jeL, jd 


= (2R) 2 ( d ~ A) - 


) 1 — d-T-d-^ b 


7T 


a r(f-A) 


r(^)r(d-A) ’ 


h(A) = f d xd v d = R3{d -A) 


3 (l+d) 




„ [s(x,y)s(y,z)s(z,x)\ A r(d)r ( 1+ ^~ A ) 3 

Working in stereographic coordinates all integrals on a sphere have the general fornf^j 


(B.l) 


I = ( 2 R) nd ~ 2 ^k<i b ki / 


d d xi 


L = i ( x + x D ai Uk<i \ x k ~ xi\ 2bkl ’ 


(B.2) 


Using rotational invariance we can always put one point to zero and multiply our integral by the 
volume of the sphere (one can always choose the point which simplifies the subsequent calculations). 
Then making iversion x / —> x//|xj | 2 we obtain (we choose x n = 0) 


I = (2 R) nd ~ 2 ^ 


k<l bkl 


2 1 ~ d 


d +1 

7r 2 


. n— 1 

n 


d d Xi 


1 


126 .- 


(B.3) 


iW) J M 

Thus we can always reduce the number of integrals by one. We want to represent each integral I 


in the Mellin-Barnes form. After this is done we can use special Mathematica packages 49,50 to 


calculate the value of the integral as a series in e 51 . All integrals I on a sphere can be represented 


diagramatically with the help of the external line and propagator (Figure [4]) . 


a 


l 


( l+x 2 


_ ]__ 

\x-y\ 2b 


x v ' x y 

Figure 4: External line and propogator for the integrals on a sphere. 

Using Feynman parameters and the Mellin-Barnes formula, P’1 we derive the following basic relations 
depicted in Figure [5} 


5 In stereographical coordintates the chordal distance is s(x, y) = 


thus I dd xV&c = I 

16 Here we mean the formula 1 

Zi + ... + Zn-i) nr=i 1 r (~Zi). 


2R\x — y | 

(l+®2)l/2( 1+?/ 2)l/2 


and the metric is ds 2 = 


(A- 1 +A' 2 + ... + A-„P T(A) ( 27 ri) 


n— 1 \rZi \r~ A — Z\ — .. 


-V-T f +l °° ... f +! °° dZt-dZn-! nr/ Xpx. 

l)n—l J— ioo J—ZOO 1 1-1=1 Z 


'r( Ad- 
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J d d xd d y ■ 

a\ b 

x y 

f d d x 

a b 

x y 

Jd d x ■ 

by! 

“< 

X 

02 7 


= r 0 (ai, a 2 , b) 


J d d x 



i +ic>o —Zl 

= j£ri) / dziT^a^zi) . 


i +*°° 

= ( 2 ^yi I dz l dz 2 dz 3 r 2 {a,b 1 ,b 2 \z l ,z 2 ,z 3 ) 

' ’ —7.rv~i 


+ioo 6 


-2l _ ~2 3 _ -2-2 

V Z 

\ ~Z\ 

' V 


= 72^6 / n dzjr 3 (a, 61,62, &3ki, _Z4 


( 27 m- . • 

v ’ —zoo 2=1 



w d—^kZ 

-z 3/ ^6 \—Z 2 


Figure 5: Basic relations for rewriting Feynman integrals on a sphere in the Mellin-Barnes form. 

where 

7 T d r(j - b)r( 0l + b - l)T(a 2 + b - l)T( ai + a 2 + b-d) 

0 ai,a2, r(|) r(ai)r(a 2 )r(ai + a 2 + 2 b - d) 


ri(a, b\z\) = 


-K d / 2 T(—z\)T(a + b — d + ^i)r (6 + 2 i)r(d — a — 2b — z{) 


r(a)r(6)r(d-a-6) 


7r d / 2 n r(-zj)r(a + 61 + b 2 - § + Yli= 1 + 


r 2 (a,bi,b 2 \zi,z 2 ,z 3 ) = 


i=l 


r(a)r(6 1 )r(b 2 )r(d-o-6i-6 2 ) 

x r(d — a — 2 &i — 26 2 — 21 — 22 — 223 ), 


-F(6i + 21 + 2 3 )r(6 2 + 22 + 23) 


71-2 n r(-2i)r(a + bi + b 2 + 63 - d + Ya =1 Zi) 


T 3 (a,bi,b 2 ,b 3 \z 1 , ...,z 6 ) = 


1=1 


-r(6i + 21 + 24 + 25) 


r(a)r(6 1 )r(6 2 )r(6 3 )r(d -a-h-b 2 -b 3 ) 

xT(b 2 + 22 + 24 + 2 6 )r(6 3 + 23 + 25 + 2 6 )r(d - a-2bi- 2 b 2 - 2 b 3 - z\- z 2 - 23 - 224 - 2z 5 — 2 z 6 ) 

(B.4) 

Using these relations one can easily write any integral on a sphere in the Mellin-Barnes form. For 
example let us consider the diagram G^\ which contributes to the free energy of </> 4 -theory at A 4 
order. This diagram is depicted on Figure [ 6 j In this graph each line represents the propagator 
Cfj,/s(x,y) d ~ 2 and we don’t include the symmetry factor, so the integral reads 


d 3) = C 8 


d d xd d yd d zd d Wy/^^y/g^y/g^ 


(B.5) 


4 ^ J s(x,y) 2 ( d 2d >s{x,z) d 2 s(x,w) d 2 s(y,z) d 2 s(y,w) d 2 s(z,w) 2 ( d 2 ) 

Working in stereographical coordinates and using rotational invarinace to set w = 0, we arrive at 
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Ill 7 . 


Figure 6: Diagram G^ . Each line represents the sphere propogator {(j)(x)<j>(y)} = C^/s(x,y) d 2 
Symmetry factor is not included in this graph. 


the following integraj^j 


Gf } = G®(2i?) 4< ^ 4_ ^ 


2 l ~ d TT— 


d d xd d zd d y 


r (^± i ) J ((1 + |x| 2 )(l + |z| 2 )(l + |y| 2 )) 4 d \x — y\ 2 ( d 2 l\x — z\ d 2 \y — z\ d 2 

(B.6) 


Diagramatically this integral is represented (up to some prefactor) in Figure [7| 


'4 — d 

iv 

= f d d xd d yd d z d — 2 / \-yr 


x 



4 - d y ^ \4 - d 

(3) 

Figure 7: Diagramatic representation for the G\ sphere integral 


Using the relations represented in Figure [5] we easily get 


=Cl(2R)^-^ 

1 r 


2 1 ~ d iT^r 

rW) 


d +1 


f lo ° d — 2 d — 2 

X . .. 3 / dz 1 dz 2 dz 3 r 2 (4:-d,d-2,——\z 1 ,Z2,z 3 )ro{4:-d-z 1 ,4:-d-Z2,^r^-z 3 ), 

J —200 ^ ^ 

(B.7) 


where the functions T 2 and Tq are given in (B.4). In practice one should try to use the relations 


in Figure [5] in a way which gives the least number of integrations over parameters Zi in the Mellin- 
Barnes representation. 


References 

[1] S. Giornbi and I. R. Klebanov, “Interpolating between a and F 1 ,” JHEP 1503 (2015) 117, 
1409.1937. 

17 We also need to make inversion to get rid of terms \x\ d ~ 2 \y\ d ~ 2 \z\ 2 ^ d ~ 2 ' 1 in the denominator. 


36 

















[2] A. B. Zamolodchikov, “Irreversibility of the Flux of the Renormalization Group in a 2D Field 
Theory,” JETP Lett. 43 (1986) 730-732. 

[3] J. L. Cardy, “Is There a c Theorem in Four-Dimensions?,” Phys.Lett. B215 (1988) 749-752. 

[4] Z. Komargodski and A. Schwinuner, “On Renormalization Group Flows in Four 
Dimensions,” J HEP 1112 (2011) 099, 1107.3987. 

[5] Z. Komargodski, “The Constraints of Conformal Symmetry on RG Flows,” 1112.4538. 

[6] C. Cordova, T. T. Dumitrescu, and X. Yin, “Higher Derivative Terms, Toroidal 
Compactihcation, and Weyl Anomalies in Six-Dimensional (2,0) Theories,” 1505.03850. 

[7] C. Cordova, T. T. Dumitrescu, and K. Intriligator, “Anomalies, Renormalization Group 
Flows, and the a-Theorem in Six-Dimensional (1,0) Theories,” 1506.03807. 

[8] I. Jack and H. Osborn, “Analogs for the c Theorem for Four-dimensional Renormalizable 
Field Theories,” Nucl.Phys. B343 (1990) 647-688. 

[9] T. Maxheld and S. Sethi, “The Conformal Anomaly of M5-Branes,” JHEP 1206 (2012) 075, 
1204.2002. 

[10] H. Elvang, D. Z. Freedman, L.-Y. Hung, M. Kiermaier, R. C. Myers, et. al, “On 
renormalization group flows and the a-theorem in 6d,” JHEP 1210 (2012) Oil, 1205.3994. 

[11] D. L. Jafferis, I. R. Klebanov, S. S. Pufu, and B. R. Safdi, “Towards the E-Theorem: Af = 2 
Field Theories on the Three- Sphere,” JHEP 06 (2011) 102, 1103.1181. 

[12] R. C. Myers and A. Sinha, “Seeing a C-Theorem with Holography,” Phys. Rev. D82 (2010) 
046006, 1006.1263. 

[13] H. Casini, M. Huerta, and R. C. Myers, “Towards a Derivation of Holographic Entanglement 
Entropy,” JHEP 05 (2011) 036, 1102.0440. 

[14] H. Casini and M. Huerta, “On the RG running of the entanglement entropy of a circle,” 

Phys.Rev. D85 (2012) 125016, 1202.5650. 

[15] H. Liu and M. Mezei, “A Refinement of entanglement entropy and the number of degrees of 
freedom,” JHEP 1304 (2013) 162, 1202.2070. 

[16] I. R. Klebanov, S. S. Pufu, and B. R. Safdi, “E-Theorem without Supersymmetry,” JHEP 
1110 (2011) 038, 1105.4598. 

[17] I. Affleck and A. W. Ludwig, “Universal noninteger ’ground state degeneracy’ in critical 
quantum systems,” Phys.Rev.Lett. 67 (1991) 161-164. 


37 








[18] T. Kawano, Y. Nakaguchi, and T. Nishioka, “Holographic Interpolation between a and F” 
JHEP 1412 (2014) 161, 1410.5973. 

[19] E. Witten, “Multitrace operators, boundary conditions, and AdS / CFT correspondence,” 

hep-th/0112258, 

[20] S. S. Gubser and I. R. Klebanov, “A Universal result on central charges in the presence of 
double trace deformations,” Nucl.Phys. B656 (2003) 23-36, hep-th/0212138, 

[21] D. E. Diaz and H. Dorn, “Partition functions and double-trace deformations in AdS/CFT,” 
JHEP 0705 (2007) 046, hep-th/0702163. 

[22] S. Giombi, I. R. Klebanov, S. S. Pufu, B. R. Safdi, and G. Tarnopolsky, “AdS Description of 
Induced Higher-Spin Gauge Theory,” JHEP 1310 (2013) 016, 1306.5242. 

[23] M. Hogervorst, S. Rychkov, and B. C. van Rees, “A Cheap Alternative to the Lattice?,” 

1409.1581. 

[24] D. Anselmi, D. Freedman, M. T. Grisaru, and A. Johansen, “Nonperturbative formulas for 
central functions of supersymmetric gauge theories,” Nucl.Phys. B526 (1998) 543-571, 

hep-th/9708042. 

[25] K. A. Intriligator and B. Wecht, “The Exact super conformal R symmetry maximizes a,” 
Nucl.Phys. B667 (2003) 183-200, hep-th/0304128. 

[26] D. Kutasov, A. Parnachev, and D. A. Sahakyan, “Central charges and U(1)(R) symmetries 
in N=1 superYang-Mills,” JHEP 0311 (2003) 013, hep-th/0308071. 

[27] J. J. Heckman, D. R. Morrison, T. Rudelius, and C. Vafa, “Geometry of 6D RG Flows,” 
1505.00009. 

[28] J. J. Heckman and T. Rudelius, “Evidence for C-theorems in 6D SCFTs,” 1506.06753, 

[29] K. G. Wilson and M. E. Fisher, “Critical exponents in 3.99 dimensions,” Phys.Rev.Lett. 28 
(1972) 240-243. 

[30] A. Hasenfratz, P. Hasenfratz, K. Jansen, J. Kuti, and Y. Shen, “The Equivalence of the top 
quark condensate and the elementary Higgs field,” Nucl.Phys. B365 (1991) 79-97. 

[31] J. Zinn-Justin, “Four fermion interaction near four-dimensions,” Nucl.Phys. B367 (1991) 
105-122. 

[32] B. Halperin, T. Lubensky, and S.-k. Ma, “First order phase transitions in superconductors 
and smectic A liquid crystals,” Phys.Rev.Lett. 32 (1974) 292-295. 


38 




[33] L. Fei, S. Giombi, and I. R. Klebanov, “Critical O(N) Models in 6 — e Dimensions,” 
Phys.Rev. D90 (2014) 025018, 1404.1094. 

[34] L. Fei, S. Giombi, I. R. Klebanov, and G. Tarnopolsky, “Three loop analysis of the critical 
O(N) models in 6- dimensions,” Phys.Rev. D91 (2015), no. 4 045011, 1411.1099. 

[35] J. Gracey, “Four loop renormalization of <f> 3 theory in six dimensions,” 1506.03357. 

[36] K. Wilson and J. B. Kogut, “The Renormalization group and the epsilon expansion,” 
Phys.Rept. 12 (1974) 75-200. 

[37] I. Drummond and G. Shore, “Conformal Anomalies for Interacting Scalar Fields in Curved 
Space-Time,” Phys.Rev. D19 (1979) 1134. 

[38] L. S. Brown and J. C. Collins, “Dimensional Renormalization of Scalar Field Theory in 
Curved Space-time,” Annals Phys. 130 (1980) 215. 

[39] S. Hathrell, “Trace Anomalies and Ac/> 4 Theory in Curved Space,” Annals Phys. 139 (1982) 
136. 

[40] I. Jack and H. Osborn, “Background Field Calculations in Curved Space-time. 1. General 
Formalism and Application to Scalar Fields,” Nucl.Phys. B234 (1984) 331. 

[41] S. El-Showk, M. F. Paulos, D. Poland, S. Rychkov, D. Simmons-Duffin, et. al., “Solving the 
3D Ising Model with the Conformal Bootstrap,” Phys.Rev. D86 (2012) 025022, 1203.6064. 

[42] S. El-Showk, M. F. Paulos, D. Poland, S. Rychkov, D. Simmons-Duffin, et. al., “Solving the 
3d Ising Model with the Conformal Bootstrap II. c-Minimization and Precise Critical 
Exponents,” 1403.4545, 

[43] H. Kleinert and V. Schulte-Frohlinde, “Critical properties of <)> 4 -theories,”. 

[44] F. Kos, D. Poland, and D. Simmons-Duffin, “Bootstrapping the O(N) Vector Models,” 
1307.6856. 

[45] D. J. Gross and A. Neveu, “Dynamical Symmetry Breaking in Asymptotically Free Field 
Theories,” Phys.Rev. DIO (1974) 3235. 

[46] M. Moshe and J. Zinn-Justin, “Quantum field theory in the large N limit: A Review,” 
Phys.Rept. 385 (2003) 69-228, hep-th/0306133. 

[47] L. Fei, S. Giombi, I. R. Klebanov, and G. Tarnopolsky, “To appear,”. 

[48] I. Drummond, “Conformally Invariant Amplitudes and Field Theory in a Space-Time of 
Constant Curvature,” Phys.Rev. D19 (1979) 1123. 


39 



[49] M. Czakon, “Automatized analytic continuation of Mellin-Barnes integrals,” 
Comput.Phys.Commun. 175 (2006) 559-571, hep-ph/0511200, 

[50] A. Smirnov and V. Smirnov, “On the Resolution of Singularities of Multiple Mellin-Barnes 
Integrals,” Eur.Phys.J. C62 (2009) 445-449, 0901.0386. 

[51] V. A. Smirnov, “Analytic tools for Feynman integrals,” Springer Tracts Mod.Phys. 250 
(2012) 1-296. 


40 




